Statistical Mechanics of the Glass Transition 
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The statistical mechanics of simple glass forming systems in 2 dimensions is worked out. The 
glass disorder is encoded via a Voronoi tessellation, and the statistical mechanics is performed 
directly in this encoding. The theory provides, without free parameters, an explanation of the glass 
transition phenomenology, including the identification of two different temperatures, T g and Tk, the 
first associated with jamming and the second associated with the appearance of a quasi-crystal at 
very low temperatures. 



Introduction: The term "glass transition" refers to 
the enormous slowing down in the dynamics of some liq- 
uids when their temperature is lowered. Despite decades 
of research, a clear explanation of this phenomenon, com- 
mon to materials as diverse as molecularglasses, metal- 
lic glasses, colloids etc., is still lacking [jlQ- The dif- 
ficulty is that the disordered molecular arrangement in 
a glass appears indistinguishable from that of the corre- 
sponding liquid, without any sign of a static correlation 
length that increases appreciably at the glass transition 
In this Letter we offer a theory of the glass transi- 
tion in simple glass forming systems, based on an encod- 
ing of the disorder that is able to flush out the pertinent 
features of the transition. In particular the disappear- 
ance of liquid-like regions, the huge increase of a typical 
scale, and eventually the existence at low temperatures 
of a quasi-crystalline state that blocks true crystalliza- 
tion 0,0, are all understood, in agreement with recent 
numerical simulations. 

The system: the glass former explored here is the 
well studied system of a two-dimensional binary mixture 
of discs interacting via a soft 1/r 12 repulsion with a "di- 
ameter" ratio of 1.4. We refer the reader to the extensive 
work done on this system |(|, showing that it is a bona 
fide glass-forming liquid meeting all the criteria of a glass 
transition. The system consists of an equimolar mixture 
of large particles with 'diameter' a-i = 1.4 and small par- 
ticles of 'diameter' o~\ = 1, respectively, but with the 
same mass m. The three pairwise additive interactions 
are given by the purely repulsive soft-core potentials 
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a,b= 1,2 , 



(1) 



where a aa = a a and a a b = {o~ a + Ob)/2. The cutoff radii 
of the interaction are set at 4.5a a b- The units of mass, 
length, time and temperature are m. <j\, r = ai^Jm/e 
and T = e/kg, respectively, with kg being Boltzmann's 
constant. A total of N = 1024 particles were enclosed in 
a square box (of area L 2 ) with periodic boundary con- 
ditions. Ref. found that for temperature T > 0.5 
the system is liquid and for lower temperatures dynami- 
cal relaxation slows down. A precise glass transition had 
not been identified in 0. 

The glass transition: To overcome the lack of sig- 
natures of the glass transition in the particle positions 



on the molecular level, we encode the state of the system 
using the Voronoi tesselation Q , where a polygon associ- 
ated with any particle contains all points closest to that 
particle than to any other particle. The edges of such 
a polygon are the perpendicular bisectors of the vectors 
joining the central particle. Such an encoding has been 
used extensively before, 0, Q, where it was also noted 
that the geometric Eulcr constraint implies that average 
coordination number is 6 at all temperatures, and local 
coordination numbers other than 6 were referred to as 
"defects" . Our encoding is richer; In our work Q it was 
discovered that a significant insight to the glass transition 
is gained by distinguishing between "liquid-like" defects 
and "glass-like" defects. We observed Q that only in the 
liquid phase there exist small particles enclosed in hep- 
tagons (or even octagons), and large particles enclosed 
in pentagons (or even squares) (cf. Fig. ^upper panel). 
In the glass phase we observe only defects of the opposite 
type, i.e. small particles in pentagons and large particles 
in heptagons, see Fig. ^ Accordingly we proposed that 
the concentration of these particular defects is a good 
indicator of the glass transition. The concentration 
q of these liquid-like defects becomes so small in 
the glass phase that we cannot distinguish it from 
zero (cf. Fig. 2] lower panel), unless the glass is put un- 
der mechanical strain, cf. Ref. Associated with this 
concentration we can define a typical scale, £: 

i = A=- (2) 



Associated with the strong decrease in ce we observe a 
huge increase in the typical scale £, in agreement with 
the tremendous slowing down of the dynamics. 

In Fig. |21 upper panel we show ci as a function of the 
temperature for a protocol of slow cooling. For tempera- 
tures larger than 0.8 the concentration follows closely an 
exponential fit, 

c e = Aexp-(AU/T) , A ps 0.094, AU ps 1.90 . (3) 

For temperatures in the range 0.3 < T < 0.8 we find an 
excellent fit to 



T g = 0.16 ±0.02 . (4) 



c e = B(T - Tg) 2 , B « 0.02, 

The quality of this fit is demonstrated in Fig. [2] lower 
panel. The fit (0} appears to identify a sharp glass tran- 
sition T g = 0.16 ± 0.02; note, however, that the fit is 
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FIG. 1: Upper panel: a Voronoi polygon construction at 
T = 3, with the seven-color code used in this paper. Small 
particles in pentagons (heptagons) are light green (dark green) 
and large particles in pentagons (heptagons) are violet (pink). 
Hexagons are in blue, octagons in black and squares in yel- 
low. Lower panel: a typical Voronoi construction in the glass 
phase at T = 0.1. Note the total disappearance of liquid-like 
defects. 

made in the range 0.3 < T < 0.8 that does not include 
T g . In fact there is no theoretical reason to expect that eg 
truly vanishes at T g , but it becomes so small that we in- 
deed do not see a single liquid-like defect in our finite-box 
simulations. Similarly, in our finite size box we cannot 
distinguish between a diverging length-scale and a scale 
much larger than L. Within the range 0.3 < T < 0.8 we 
fit an apparently divergent length 

£~(T-T g )-', u=l. (5) 

Finally, it was discovered in || that at sufficiently low 
temperatures, there exists a quasi-crystalline phase, see. 
Fig. |21 This phase arises while heating, starting from 
a square lattice with interchanging small and large par- 
ticles. The quasi-crystalline phase appears stable in our 
slow heating simulations as long as the temperature is 
sufficiently small. We propose to identify the quasi- 



0.07 
0.06 
0.05 
0.04 

cT 

0.03 
0.02 
0.01 



0.5 1 1.5 2 2.5 3 

T 




FIG. 2: Concentration of defects under slow cooling. Up- 
per panel: the concentration of the liquid-like defects: large 
particles in pentagons (red dots) and small particles in hep- 
tagons (blue dots). Lower panel: the fit of the concentration 
of liquid-like defects as a function of temperature according 
to Eq. JIJ. 




FIG. 3: The 'ideal' glass state. 

crystalline phase as the ideal glass state, and the tem- 
perature where it melts as the Kauzmann temperature 
Tfc 0. Upon heating above Tk the quasi-crystal melts 
in favor of a generic glassy state in which the quasi- 
crystalline phase is invaded by patches of hexagons, see 
Fig. flower panel. To understand these findings we need 
to understand why the liquid-like defects disappear first 
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upon cooling, and what is the thermodynamic reason for 
the non-existence of hexagons at sufficiently low temper- 
atures. 

Statistical Mechanics: here we present the statisti- 
cal mechanics of this system, rationalizing and explaining 
the features of the transition presented above. The main 
point to stress is that even in the glass phase, when par- 
ticles are practically jammed and precluded from large 
excursions, the Voronoi tesselation is still highly mobile. 
Minute changes in particle positions can lead, via Tl 
foam-like processes [l(| to transitions between cell-types. 
Therefore, if one can convince oneself that the cell types, 
including the presence of a small or large particle inside 
them, are proper 'species', one can construct statistical 
mechanics right on the space of these cells. We define the 
energy of each cell type as the average (over all cells of the 

E ( \ 12 x 

same type) of the potential energy ej = Q^fc=i e ( 7^ ) ) 
where Ei is the number of edges associated with that cell 
type, and being the distance to the particle in the 
adjacent Voronoi cell. In Fig. 0] we present the values 
of these energies measured as a function of the tempera- 
ture, following a protocol of slow cooling. There are 10 
different cell types in this system, (large paricle or small 
particle in squares, pentagons, hexagons, heptagons and 
octagons), but octagons and squares are not shown since 
they disappear much before the glass transition. The 
lesson drawn from this graph is that the different cell 
types have clearly split energies throughout the 
interesting temperature range, and that these en- 
ergies are only weakly dependent on the temper- 
ature. Within the temperature range of interest we can 
focus on the six types of cells; denote by {iVj}| =1 the 
number of cells of each type, with number of edges Ei, 
ordering them by the energy e 2 ; from i = 1 being the 
highest one (large particle in a pentagon) to i = 6 being 
the small particle in an heptagon. Additional important 
properties of the cell types are their areas Qi and their 
shapes; the first affects the enthalpy term and both af- 
fect the entropy when we count the number of possible 
tilings of the plane. With this in mind we can construct 
the statistical mechanics of this system by considering 
the free energy G = U +pV — TS. The value of U is then 

simply J2i=i N i £ i- Tne P v term is simply pJ2i=i 
Lastly, we need to estimate the entropy term. In princi- 
ple this should be computed from the number of possible 
complete tilings of the area by iVj cells of each type with 
its given area and shape, subject to the Eulcr constraint 
Y^i—i NiEi = 6N. This is a formidable problem. A use- 
ful estimate can be made by considering the area only, 
and filling space starting with the largest objects, then 
the next largest etc. until the smallest is fit in. Denot- 
ing the possible number of boxes to fit the largest cells 
by Afi = V/fti, then the number of boxes available for 
the second largest cell by A2 = (V— Nitti)/^ etc., the 
number of possible configurations W is 
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• 1=1, large in pentagons 
A 1=2, large in hexagons 

• i=3, large in heptagons 

• 1=4, small in pentagons 
A 1=5, small in hexagons 

1=6, small in heptagons 
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FIG. 4: The average energies of the voronoi cells as a function 
of the temperature as measured in the simulations. 

Denoting xi = Ni/Afi we compute directly Xi = 
CjQj/ Ylj—i CjClj where c,; is the number concentration 
of each defect. We can now compute S = fee In W and 
write G together with a Lagrange multiplier for the Eulcr 
constraint, 

6 6 6 

G = J2 N id + P NiQi + X J2 NiEi 

i—1 i—1 i—1 

6 

+Tj2^k[x k ^x k + (l-x k )lii(l-x k )] . (7) 
fe=i 

The chemical potential fii = dG/dNi is then 

m = e t +pSli + T[ln Xl +J2cT ln ( X ~ Xk ^ + XE * ■ ( 8 ) 

fe=l k 

We now recognize that in equilibrium there exist only 
two independent values of /i;, one for the small particles 
Hs and one for the large particles hl, and we have 9 
unknowns - six values of c^, 2 values of fj. and one La- 
grange multiplier A. This is precisely balanced by the 
6 equations (JSJ), the Euler constraint, and the two con- 
straints J2i=i c i = 12k=4 c k = 1/2- These equations can 
be solved numerically using the precise values of f2j(T) 
and 6i(T) as measured in the simulation. The approxi- 
mate calculation of the entropy however does not warrant 
such a detailed calculation. In reality, calculating the av- 
erage areas of the cell types in the numerical simulations, 
we discover that to an excellent approximation these fall 
in two classes, smaller cells of area £ls when small parti- 
cles are in them, and larger cells of area £1^, when large 
particles are enclosed. These areas again are only weakly 
dependent on the temperature. Then the whole system 
of equations simplifies to two analytically tractable sets 
of equations 

pi L = ^ + T In c, t + \E, t , { * = 1 , 2, 3} , 

ft s = et+ThiCi + XEi , {z = 4,5,6}, (9) 
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FIG. 5: The concentrations of Votonoi cells as predicted by 
the statistical mechanics. Note the close-to-degeneracy of 
pairs of cells: the upper pair are the glass-like defects, the 
middle pair are the two hexagons and the lower pair the liquid- 
like defects. Observe the existence of two temperatures, in- 
terpreted as T g and Tk . Lower panel: the concentration of the 
liquid-like defect (small in heptagon) as a function of temper- 
ature, together with the quadratic fit of the type Eq. 

together with the above mentioned three constraints. In 
p, we have absorbed terms that added to fj, in this special 
case. In Fig. [S]we show the solutions of these equations 
when we use values of taken from the Fig. 0] at T = 
. The results are in excellent agreement with the sim- 
ulations. There is a clear degeneracy in the temperature 
dependence of the two types of liquid like defects, the 
two hexagons and the pair of glass-like defects. At low 
temperatures we only have the glass-like defects who can- 
not crystalizc, but can form the quasi-crystalline phase. 



Upon warming up, at TJ-, a sizable number of hexagons 
appears to form the generic glassy state. Yet later, at T g , 
a sizable number of liquid-like defects brings the system 
to a liquid state. The curves associated with the tem- 
perature dependence of the concentration of the latter 
defects can be very well fit by a quadratic fit like Eq. 
(J3J, with a temperature T g « 0.07. In fact, the solution 
of Eqs. © shows that the concentrations do not quite 
vanish but become exponentially small, as stated above 
and in ||. Note the numerical value of T g differs from 
the simulation, as is expected from a theory that does not 
take careful account of the correlations between different 
cells. The actual values of T g and Tk can be understood 
from this model. Denote by ce, ch and cg the concentra- 
tions of the liquid-like, hexagons and glass-like defects, 
and by = e\ + eg ~ 12.48 as the energy associated with 
the liquid- like defects, by eh = £2 + £5 ~ 11.94 the en- 
ergy of the hexagons, and eg = £3 + 64 ~ 11.76 the energy 
of the glass-like defects. The theory predicts that ratios 
ci/ch and ch/cq are of the order of exp[— (eg — e#)/T] 
and exp[— (e# — €q)/T] respectively. As an estimate of 
T g and Tk take these ratios to be, say, of the order of 
1% ~ exp(— 5) and observe that such ratios are obtained 
for T = T g w 0.11 and T = T K w 0.04. It is important to 
notice that e#— eg could be positive rather than negative, 
and then the system would crystallize on a hexagonal lat- 
tice! 

Conclusions: It is quite remarkable that one can un- 
derstand the properties of a system of strongly interact- 
ing particles in an essentially jammed state using the 
statistical mechanics of a mobile set of essentially non- 
interacting quasi-particles. Without any free parameter 
we have offered a semi-quantitative understanding of the 
features of the glass transition in this system, as well as 
the existence of two distinct temperatures T g and Tk. 
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